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Abstract. 

During the year 2000 the CERES experiment, upgraded with a radial drift 
TPC, took a large data sample of Pb on Au collisions at 158 AGeV triggered on the 
8% most central collisions. A very detailed calibration of the radial drift TPC was 
since completed. Preliminary results on e+e^ pairs and if> mesons reconstructed 
(~^ ' in the K^K~ channel are presented. 
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1. Introduction 



^ ' QCD calculations predict a transition from ordinary hadronic matter into a plasma of 

deconfined quarks and gluons at high energy density. Dileptons, which have negligible 
^^ ' final state interactions, represent a very suitable probe for the study of this new 

\Q '. state of matter. The CERES/NA45 experiment at the CERN SPS has made major 

^^ ' contributions to the measurement of low-mass e"'"e~ pairs in ultra-relativistic heavy- 

^r , ion collisions. CERES has measured an enhanced dilepton production in the invariant 

mass region rrii.+ e'- > 0.2 GeV/c^ in S+Au at 200 AGeV T and in Pb+Au at 



o 



o 

^ , 158 AGeV |2j and at 40 AGeV PI compared to the contribution from known hadronic 

^ ' sources. The enhancement is absent in p-induced reactions 0]. Pion annihilation has 

been taken into account as an additional mechanism for e"'"e~ production but the 
r-{ ' experimental spectra cannot be explained without introducing medium modifications 

^ , of vector mesons, particularly of the p. The question of possible modifications of other 

vector mesons, specially of the w and the 0, and the role of chiral symmetry restoration 
remained however open. 

In order to further investigate the enhancement, the CERES spectrometer was 
upgraded during 1998 by the addition of a Time Projection Chamber (TPC) with 
radial electric drift field [51 [Q which improves the mass resolution and the electron 
identification. During the year 2000 CERES took a large data sample consisting of 
30T0^ and 3T0^ events of Pb on Au collisions at 158 AGeV triggered on the 8% and 
20% most central collisions, respectively. The challenge was to calibrate the highly 
complex radial drift TPC to bring the momentum resolution to the design limit. 

2. Experimental Setup 

The CERES experiment (Fig. ^ is optimized to measure low mass electron pairs 
close to mid-rapidity (2.1< rj <2.6) with full azimuthal coverage. A vertex telescope, 
composed of two Silicon Drift Detectors (SDD) positioned at 10 cm and 13.8 cm 

t For the full CERES Collaboration author list and acknowledgments, see Appendix "Collaborations" 
of this volume. 
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Figure 1. Cross section through the upgraded setup of the CERES spectrometer 

downstream of a segmented Au target, provides a precise vertex reconstruction, angle 
measurement for charged particles and rejection of close pairs from 7 conversions and 
tt" Dalitz decays. Two Ring Imaging CHerenkov (RICH) detectors, operated at a high 
threshold (74^=82), are used for electron identification in a large hadronic background. 
The new radial-drift TPC, positioned downstream of the original spectrometer, has 
an active length of 2 m and a diameter of 2.6 m. A gas mixture of Ne (80%) and 
CO2 (20%) is used. It is operated inside a magnetic field {B, indicated by the dashed 
field lines inside the TPC in Fig. QJ with a maximal radial component of 0.5 T and 
provides up to 20 space points for each charged particle track. This is sufficient for the 
momentum determination and for additional electron identification via d£'/da; in the 
TPC. In the configuration with the TPC, the magnetic field between the two RICH 
detectors is switched off and there is thus no deflection between them, allowing to use 
them in a combined mode resulting in an increased electron efficiency. Moreover, the 
TPC opens new possibilities to study hadronic observables |51l^ 1101 1111 112| . 

3. Calibration of the TPC 



One of the main tasks of the TPC calibration was to understand the electric field 
{E) which is dominantly radial with Er ~ 1/r but has a small azimuthal component 
due to the polygonal shape of the TPC and small longitudinal components at the end 
caps. A custom program based on the relaxation method was developed for calculating 
the electric field including the proper knowledge of the field cage resistors, the field 
distortions caused by slightly displaced chambers and the leakage of the amplification 
field through the cathode wire plane. 

The magnetic field in the TPC is generated by two coils running with opposite 
currents (Fig^l. The two main components Bz and B^ change with r and z. A 
field map was measured before the installation of the TPC in the experimental area. 
Compared to the nominal magnetic field calculated with the POISSON program the 
measured field deviates from azimuthal symmetry by a few %. These deviations are 
now included as corrections to the nominal field map. 
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The electron mobility as a function of E was determined from laser events 
containing tracks at known positions for a given set of gas parameters. Relative 
run-by-run variations caused by changes in composition, temperature, and pressure 
(monitored by the slow control system of the experiment) are calculated with the 
MAGBOLTZ program ^31- The transformation from pad, time, and plane to the 
laboratory coordinate system (x, y, and z) is done using a fourth-order Runge-Kutta 
method. The drift trajectory is calculated starting at the cathode plane using in 
each point the drift velocity vector, vd — /^ i+( \b\)'^ ( '^ 1^^ x B) + fi'^{E ■ B)B) 

where the components parallel to the electric field and parallel to {E x B) have been 
modified to account for the difference observed between using vd from the above 
equation and the Monte Carlo drift option in the MAGBOLTZ program J^!- To 
fully account for the Lorentz angle a correction is extracted from the symmetry of the 
1/p distribution of identified pions under the assumption that very high momentum 
particles have vanishing curvature. After all these corrections the local distortions are 
small compared to the resolution. 

The momentum of tracks from the target is determined from the combination 
(resolution-weighted) of a 3-parameter fit that takes into account multiple scattering 
before the TPC (dominant for low-momentum particles) and a 2-parameter fit that 
assumes that tracks come from the vertex (dominant for high-momentum particles). 
The momentum resolution depends on the number of hits on a track and on the single- 
hit position resolution. New hit finding and track finding algorithms were developed 
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Figure 2. (Left) Single-hit position resolution in azimuthal direction (rA(/>) 
as function of the plane in the TPC (z) for events triggered on the 8% (circles) 
and 20% (triangles) most central collisions obtained from hits on tracks with a 
momentum p > 1 GeV/c. (Middle) Momentum resolution from a simulation with 
single hit position resolution as in the data for the 3 parameter fit (triangles), for 
the 2 parameter fit (squares) and for the combined fit (circles). (Right) Invariant 
mass spectra of p7r~ (pvr+) pairs after background subtraction for 2.0 < VA < 2.4 
and pj > 1.5 GeV/c. After fitting the peaks of A (circles) and A (squares) with 
a Gaussian function widths of 6.94±0.02 MeV/c^ and 6.85±0.13 MeV/c^ are 
obtained. The measured A/A ratio is 0.105±0.002. 



to better resolve close hits and to improve tracking efficiency at very low momenta. 
Tracks above 1 GeV/c momentum are found with an efficiency of 90% due to non 
instrumented areas and have an average number of hits of 18.7. Down to a momentum 
of 0.3 GeV/c the efficiency is still 80% and the average number of hits is 17.2. The 
single-hit position resolution is obtained comparing the reconstructed hit positions 
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with the ideal hit positions given by the fitted trajectory of a track. The width of 
the residual distribution gives the spatial resolution. The resolution in azimuthal 
direction scaled with the radius r of the hit (urA.^) is shown in Fig.|21 It deteriorates 
slightly with increasing hit multiplicity in the TPC. The momentum resolution has 
also been evaluated using a microscopic drift Monte Carlo simulation with single- 
point resolution as observed in the data (Fig. El . It translates into a mass resolution 
of about 4% at the meson mass reconstructed in the e^e~ channel. The momentum 
resolution was checked in the data by looking at the width of the reconstructed A and 
K^ in the invariant-mass spectrum (Fig. O • The detailed calibration improves the 
mass resolution by at least a factor 1.7 compared to the old cahbration ^^. 
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Figure 3. (Left) dE/dx resolution of electrons as a function of number of hits in 
the track after the new calibration (squares) compared to the parameterization 
of Allison/Cobb (line). (Right) dB/dx signal in the TPC versus momentum for 
electrons as selected with the RICH detector. The inclined cut against pions 
(black line) show the additional rejection power from the TPC. 

The pad-to-pad gain calibration is obtained from the total hit amplitude in the 
pad of maximum amplitude. A correction of the undershoot after each pulse obtained 
from non-zero suppressed data was included [15!. An electron attachment correction 
is determined by analysing the dependence of the hit amplitude on the drift length 
and taking into account the different particle abundancies for each polar angle 0. The 
dii'/dx of a track is calculated as the truncated mean of the individual hit amplitudes. 
For the maximum number of hits on tracks the d-E/dx resolution approaches the 
parameterization by Allison and Cobb (Fig.|3J), which implies that most of the electrons 
in this analysis are measured with di?/dx resolution of better than 10%. 

4. Electron analysis 



The main difficulties of the electron analysis are the low probability of electromagnetic 
decays and the large amount of combinatorial background from 7 conversions and 
Dalitz decays. To select electrons among all charged hadrons the RICH detectors and 
the diJ/dx signal in the TPC are used. The analysis chain provides track segments in 
the different detectors. The two RICH detectors are used in a combined mode in this 
analysis. Rings with asymptotic radius and a (double Hough) amplitude above a given 
threshold for the full ring acceptance (slowly decreasing threshold towards the edges) 
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are taken as RICH track segnients. With these cuts a pion rejection of about 0.999 
is achieved for an electron efficiency of about 0.55. QuaUty cuts (matching between 
SDDl/2 of 2 mrad and at least 12 hits in the TPC track segments) are applied. 

Global tracks are constructed by combining track segments in SDD, TPC and 
RICH with a 2a momentum-dependent matching window. The electron d_E/dx signal 
in the TPC is independent of momentum (Fig. |3I). A primary electron selection is 
done based on the d-E/dx resolution. Finally a more restricted cut to further suppress 
pions is applied (Fig.OJ. 

Very good electron identification is not enough. Electron pairs from 7 
conversions and i:^ Dalitz decays need to be identified by topology and removed 
from the sample to reduce combinatorial background. They are characterized by 
their small opening angle and low momentum. As the detectors have a finite 
two-track resolution the most effective way of rejecting conversions in the target 
and close Dalitz pairs is by rejecting tracks with at least twice the minimum 
ionizing energy loss signal in both SDD's where all hit amplitudes have been re- 
summed within 7 mrad around the each SDD track segment. Low- amplitude tails 
in the SDD di?/dx distribution are also removed. Late conversions (mostly in 
SDDl) are removed by an upper cut in each RICH ring amplitude and by a 
cut in the distance to another TPC track of opposite sign and electron d-E/dx. 
To reduce 7r° Dalitz decay contributions 
only electron tracks with no opposite 
charge electron track within 35 mrad are 
taken for analysis. Finally, only electron 
tracks in the geometrical acceptance 
0.14 rad < 6 < 0.243 rad and with a 
Pt > 0.2 GeV/c are selected. 

The invariant-mass distributions of 
unlike-sign pairs and of like-sign pairs 
that is a measure of the combinato- 
rial background after full rejection are 
shown in Fig. 21 In order to reduce 
the statistical errors an unlike-sign com- 
binatorial background using the mixed- 
event technique has also been evalu- 
ated. For masses below 0.2 GeV/c^, the 
mixed-event background deviates from 
the like-sign background due to ring 
reconstruction effects; the mixed-event 
background is normalized to the like- 
sign pair background in the mass region 




iTiee (GeV/c ) 



Figure 4. Invariant-mass dis- 
tribution of unlike-sign pairs 
(full symbols), like-sign pairs 
(open symbols), and mixed- 
event background (line) nor- 
malized to like-sign pairs back- 
ground after full rejection. 



0.2 GeV/c^ <me+e- <0.6 GeV/c^ The 
physics signal is obtained by subtract- 
ing the like-sign pairs or the mixed event 
background from the unlike-sign pairs as 

Ng+g- — (7Ve+e+ -I- A^e-e- )• This analysis is based on about 18-10^ events corresponding 
aproximately to 80% of the events where both RICH detectors were operational. 
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5. Low-mass e+e results 

The preliminary e^e~ invariant-mass spectrum normalized to the tt*^ Dalitz peak 
of the hadronic decay cocktail is shown in Fig. The hadron decay cocktail [T7] 
has been folded with the experimental momentum resolution and energy loss due 
to bremsstrahlung. Acceptance, opening-angle, and transverse-momentum cuts are 
applied. An excess of pairs for m^+g- > 0.2 GeV/c^ is visible. Using the like- 
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Figure 5. Invariant-mass spectrum of e+e~-pairs normalized to the expectation 
from the hadron decay cocktail (thick solid line) in the n^ Dalitz peak using 
the like-sign background (left) or the mixed-event background (right). The 
expectations from model calculations assuming the vacuum p spectral function 
(dashed), a dropping mass (dotted) or a in medium spread p width (dashed- 
dotted) are also shown. 

sign background, the number of pairs in the Dalitz region (mg+g- < 0.2 GeV/c^) 
is 2826 ± 114 with a signal to background ratio (S/B) of 1/1.78. Using the mixed- 
event background, the number of open pairs (m^+g- > 0.2 GeV/c^) is 2023 ± 182 
with a S/B of 1/14.6. The enhancement factor for 0.2 GeV/c2< mg+g-<l.l GeV/c^ 
compared to the hadron decay cocktail is 3.1 ± 0.3 (stat). Further optimization of 
cuts based on Monte Carlo simulations is in progress and is expected to improve the 
statistical errors. 

The experimental results are compared to theoretical models based on hadronic 
decays and 7r"''7r~ annihilation. The p-propagator is treated ^S] in 3 ways: vacuum 
/9, modifications following Brown- Rho scaling |19| . and modifications via p-hadron 
interactions [201 ■ The data clearly disfavour the vacuum p. In the region between 
the to and the the data seem to favour the many-body approach over Brown-Rho 
scaling. We have to await the final results of the analysis before a stronger statement 
can be made. 



6. The (j) meson puzzle 

One of the open issues at the CERN SPS is the (f> puzzle. Two experiments, one 
measuring the (j) meson in the /x"*" fi~ decay channel (NA50), and the other (NA49), 
measuring it in the K^K~ channel, obtain yields that differ by factors between 2 and 
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4 in the common m^ range |21| . Further, the rrit spectra exhibit a different inverse 
slope parameter [HI, 305±15 MeV in NA49 and 2f8 ± 6 MeV in NA50, fitted in their 
rrit acceptance regions. 
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Figure 6. Invariant-mass spectrum of e+e~-pairs together witli tlie expectations 
from model calculations assuming in-medium spread p width plus the thermal 
yield (thin solid) or the thermal yield reduced to 50% (dashed) are shown. The left 
plot corresponds to the standard set of cuts and the right one has an additional 
electron pair pt cut of 0.3 GeV/c. 



We have compared our preUminary mass spectrum with a calculated spectrum 
based on a (/)-yield as predicted in a thermal model [221 which also describes the NA50 
yield (left of Fig.|HJl. The lines also contain the p-contribution. Within a window of 
±lcr around the </> peak the ratio of data to simulation is 1.1 ± 0.3. With an additional 
electron pair pt cut of 300 MeV/c (right of Fig. O the ratio is 1.2 ± 0.3. Obviously 
this is in very good agreement with our data. Reducing the (/)-yield to 50% of the 
above value (dashed hues of Fig. EJ) leads to a ratio of 1.9 ± 0.5 and 2.1 ± 0.6. The 
reduced yield that would agree with the NA49 observation for </> -^ K^K~ seems in 
disagreement with our preliminary data. The evaluation of absolute efficiencies and 
systematic errors are needed to draw a definitive conclusion. 

The upgraded spectrometer offers the unique possibility to study simultaneously 
the 4) meson in the K^K~ decay channel. To study this channel, we assign to 
all charged particles the kaon mass (no particle identification). TPC tracks in 
the geometrical acceptance 0.14 rad < 6 < 0.24 rad with at least 12 hits and 
with a transverse momentum pt > 0.25 GeV/c that are matched to a SDD track 
within 2.5(7 momentum-dependent matching window are used. The combinatorial 
background is calculated using the mixed-event technique. An opening-angle cut of 
0.015 rad < 0i2 < 0.17 rad is applied to improve the signal-to-background ratio. A 
cut qt < 0.19 — 0.4 • a"* — 0.4 • a^, where a is the Podolanski-Armenteros parameter 
a — (pj — P~[)/{pt, + Pl) ^^'^ Pl (it) i^ the decay track momentum component 
parallel (transverse) to the momentum in the laboratory frame, is also applied to 
avoid contamination from other particles species. The signal is studied for pf > 1 
GeV/c and rapidity 2.0 < y"^ < 2.4. Acceptance and pair reconstruction efficiency 
corrections are obtained from GEANT simulations where (j) decays with a distribution 
in transverse momentum dN/dpt = C ■ pt ■ exp{—mt/T) and a Gaussian rapidity 
distribution with a standard deviation of 1.2 are embedded into the background of 
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real events. The preliminary invariant mass spectrum after background subtraction 
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Figure 7. Preliminary invariant-mass spectrum of K^K pairs after background 
subtraction for 2.2 < y^ < 2.4 and 1.5 GeV/c < pf < 1.75 GeV/c . A peak 
corresponding to the meson is observed (left). Transverse momentum spectrum 
of the <j) meson corrected for acceptance and efficiency. Also shown are fits with 
a mentioned distribution. For fits see text. 

and the transverse momentum spectra corrected for efficiency and acceptance are 
presented in Fig. [7| The values of the inverse slope parameter obtained after fitting 
with the function given above are T=218±12 MeV and T=246±22 MeV for the 
rapidity intervals 2.0<y<2.2 and 2.2<y<2.4, respectively. Evaluation of systematic 
errors and of the exact centrality necessary to compare to other experiments are in 
progress. 
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